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Abstract. We develop the max-plus finite element method to solve finite 
horizon deterministic optimal control problems. This method, that we intro- 
duced in a previous work, relies on a max-plus variational formulation, and 
exploits the properties of projectors on max-plus semimodules. We prove here 
a convergence result, in arbitrary dimension, showing that for a subclass of 
problems, the error estimate is of order S + Ax{5)~^ , where 5 and Ax are the 
time and space steps respectively. We also show how the max-plus analogues 
of the mass and stiffness matrices can be computed by convex optimization, 
even when the global problem is non convex. We illustrate the method by 
numerical examples in dimension 2. 



1. Introduction 
We consider the optimal control problem: 

(la) maximize / £{x{s),u{s)) ds + (f>{x{T)) 

Jo 

over the set of trajectories {x{-),u{-)) satisfying 

x{s) = f{x{s),u{s)), x{0)=x, 
^ ' x{s) e X, u{s) e U, 

for all < s < T. Here, the state space X is a subset of R", the set of control 
values [/ is a subset of R"*, the horizon T > and the initial condition x (z X 
are given, we assume that the map u{-) is measurable, and that the map a;(-) is 
absolutely continuous. We also assume that the instantaneous reward or Lagrangian 
£ : X X U ^ R, and the dynamics f : X x U — > R", are sufficiently regular maps, 
and that the terminal reward is a map X — > R U {— cx)}. 

We are interested in the numerical computation of the value function v which 
associates to any {x,t) £ X x [0,T] the supremum v{x,t) of £{x{s),u{s)) ds + 
(j){x{t)), under the constraint Ijlbp . for Q < s < t. It is known that, under certain 
regularity assumptions, v is solution of the Hamilton- Jacobi equation 

3v 3v 

(2a) __+ij(a;,_) = o, (x,t)eXx (0,T] , 

with initial condition: 

(2b) v{x,Q) = (j){x), xeX , 
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where H{x,p) = sup^^^u £{x,u) + p ■ f{x,u) is the Hamiltonian of the problem (see 
for instance |Lio82l IBar94j ). The evolution semigroup S** of |5J), or Lax-Oleinik 
semigroup, associates to any map tj) the function := f(-,t), where v is the value 
function of the optimal control problem 

Maslov observed in |Mas73) that the evolution semigroup S** is max-plus linear 
(see also jMS92l IKM97j V Recall that the max-plus semiring, Kmax, is the set 
R U {—(X)}, equipped with the addition a © & = max(a,6) and the multiplication 
a ® h = a + h. By max-plus linearity, we mean that for all maps /, g from X to 
Rmax, and for all A G Mmax, we have 



where f®g denotes the map x ^ f{x)®g{x), and A/ denotes the map x ^ X®f{x). 
Linear operators over max-plus type semirings have been widely studied, see for 
instan ce IHnTOl IMS92I [BCODQ"^ IKMH7I IC^MOU . see also |FataF) . 

In |FMOO| ■ Fleming and McEneaney introduced a first discretization method 
exploiting the max-plus linearity of the semigroup . 

In |AGL04j . we introduced a new max-plus based discretization method, in- 
spired by the classical finite element method. The max-plus finite element method 
of | AGL04j approximates the evolution semigroup 5* by means of a nonlinear dis- 
crete semigroup, which can be interpreted as the dynamic programming operator 
of a deterministic zero-sum two players game, with finite action and state spaces 
(unlike the method of Fleming and McEneaney which leads to a discrete optimal 
control problem). The state of the game corresponds to the set of finite elements. 
To each test function corresponds one possible action of the first player, and to 
each finite element corresponds one possible action of the second player. This 
discretization, which can be interpreted geometrically in terms of projections on 
semimodules, is similar to the classical Petrov-Galerkin finite element method. 

The computation of the instantaneous payments of the game requires the eval- 
uation of the max-plus scalar product {z \ S^w) for each finite element w and 
each test function z, where 5 is the time discretization step. In some special cases, 
{z I S^w) can be computed analytically. In general, we need to approximate this 
scalar product, for each finite element w and test function z. In )A(TL04j . we used 
the simplest approximation S^w = w -\- SH{-, ^), already considered in |MH99j . 
This requires regularity assumptions on w (or alternatively, on z, if one uses the 
dual semigroup |AGL04| '). In this paper, we rather use a direct method, which 
allows us to approximate {z \ S^w) by the value of an optimization problem in 
finite dimension. We show that, under reasonable assumptions on £, /, z and w, 
this approximation leads to a concave optimization problem. We also give an error 
estimate of order S -\- 

The paper is organised as follows. In Section |21 we recall some basic tools 
and notions: residuation, semimodules and projection. In Section 13 we recall the 
formulation of the max-plus finite element method. The contents of Sections [21 
and O are essentially taken from |AGL04j : we need to recall them to state our 
results. In Sectional we discuss the approximation of the scalar product {z \ S^w). 
In Section we give the main convergence theorem. Finally, in Section El we 
illustrate the method by numerical examples in dimension 2. 



S'{f®g) 
S'iXf) 



S'f © S*g , 
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2. Preliminaries on residuation and projections over semimodules 

In this section we recall some classical residuation results (see for example 
|B J 72) . |BCOQ92 ), and their application to linear maps on idempotent semimod- 
|C(JQ04| ). We also review some results of |CUQ96| |CCiQ04| 



ules (see |LMS()1 



concerning projectors over semimodules. 

2.1. Residuation, semimodules, and linear maps. If {S,<) and (T, <) are 

(partially) ordered sets, we say that a map / : S* — > T is monotone if s < s' =J> 
f{s) < /(s')- We say that / is residuated if there exists a map f^:T^S such 
that 

f(s) <t ^ s< fit) . 

The map / is residuated if, and only if, for all t E T, {s E S \ f{s) < t} has a 
maximum element in S. Then, 

fit) = max{s G S \ f{s) <t}, Vt € T . 

If a set /C is a monoid for a commutative idempotent law © (idempotent means that 
a ® a = a) , the natural order on JC is defined by a < 6 <^=J> a 6 = 5. We say that 
/C is complete as a naturally ordered set if any subset of JC has a least upper bound 
for the natural order. If (/C, ®, (8") is an idempotent semiring, i.e., a semiring whose 
addition is idempotent, we say that the semiring /C is complete if it is complete as 
a naturally ordered set, and if the left and right multiplications, L^, i?^ : /C ^ /C, 
L^(x) = ax, R'^ix) = xa, arc residuated. 

The max-plus semiring, Mmax; is an idempotent semiring. It is not complete, but 
it can be embedded in the complete idempotent semiring Rmax obtained by adjoin- 
ing +00 to Mmax, with the convention that —00 is absorbing for the multiplication 
a®h = a + h. The map x 1— > —x from M to itself yields an isomorphism from Mmax 
to the complete idempotent semiring Mminj obtained by replacing max by min and 
by exchanging the roles of +00 and — cx) in the definition of Mmax- 

Semimodules over semirings are defined like modules over rings, mutatis mutan- 
dis, see |LMS01I|CGQ04| . When /C is a complete idempotent semiring, we say that 
a (right) /C-semimodule X is complete if it is complete as a naturally ordered set, 
and if, for &\\ u E X and A S /C, the right and left multiplications, i?^ : X ^ X , 
V i-^ v\ and Ljf : /C — > A", /i i— > w/i, are residuated. In a complete semimodule A", 
we define, for all u,v E X, 

u\v =^ (L^)Hv) = max{A £ JC \ u\ < v} . 

We shall use semimodules of functions: when X is a set and {JC, ®, (81) is a complete 
idempotent semiring, the set of functions JC'^ is a complete /C-semimodule for the 
componentwise addition (u,v) 1-^ u (B v (defined by {u © v){x) = u{x) © v(x)), and 
the componentwise multiplication (A, u) 1— > u\ (defined by {uX){x) = u{x) (8 A). 

If JC is an idempotent semiring, and if X and y are /C-semimodules, we say that 
a map A : X ^ y is, additive if for all u,v £ X, A{u © w) = A{u) © A{v) and 
that A is homogeneous if for all u S A" and X G JC, A{uX) = A{u)X. We say that 
A is linear, or is a linear operator, if it is additive and homogeneous. Then, as in 
classical algebra, we use the notation Au instead of A{u). When A is residuated 
and V € y, wc use the notation A\v or A^v instead of A^{v). 

If X and Y are two sets, {JC, ©, (g)) is a complete idempotent semiring, and 
a G JC^^"^ , we construct the linear operator A from JC^ to JC-^ which associates to 
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any u G JC^ the function Au G K,^ such that Au{x) = \/ yeY a{x, y)(E)u{y), where V 
denotes the supremum for the natural order. We say that A is the kernel operator 
with kernel or matrix a. We shall often use the same notation A for the operator 
and the kernel. As is well known (see for instance |BCOQ92| ), the kernel operator 
A is residuated; and 

iA\v)iy) = A A{x,y)\vix), 

where A denotes the infimum for the natural order. In particular, when K, = Kmax, 
we have 

(3) {A\v){y) = \{-A{x, y) + v{x)) = [-A*{-v)]{y) 

x^X 

where A* denotes the transposed operator JC^ —> /C^, which is associated to the 
kernel A*(y, x) = A{x, y). (In (j2Jl, we use the convention that +00 is absorbing for 
addition.) 

2.2. Projectors on semimodules. Let V denote a complete suhsemimodule of a 
complete semimodule X over a complete idempotent semiring /C, i.e., a subset of 
X that is stable by arbitrary sups and by the action of scalars. We call canonical 
projector on V the map 

(4) Pv:X^X, u Pv(w) = max{w e V I w < It}. 

Let W denote a generating family of a complete suhsemimodule V, which means 
that any element w G V can be written as v = \/ {wX^ \ w G W}, for some X^j G /C. 
It is known that 

Pviu) = V w(w\u) 

(see for instance |CGQ04| ). li B : U ^ X is a residuated linear operator, then the 
image im_B of i? is a complete subsemimodulc of X, and 

(5) Pn-nB^BoBK 

The max-plus finite element methods relies on the notion of projection on an im- 
age, parallel to a kernel, which was introduced by Cohen, the second author, and 
Quadrat, in |CGQ96| . The following theorem, of which Proposition [21 below is an 
immediate corollary, is a variation on the results of |CGQ96[ Section 6] . 

Theorem 1 (Projection on an image parallel to a kernel). Let B : U ^ X and 

C : X y be two residuated linear operators. Let 11^ = B o [C o B)^ o C . We 
have = Eg o n'^, where Ug = B o B^ and 11*^ = o C. Moreover, Ilg is a 
projector ((11^)^ = H^), and for all x G X: 

Ub{x) = max{y G imS | Cy < Cx}. 

The results of |CGQ96| characterize the existence and uniqueness, for all x € X , 
oi y G iuiB such that Cy = Cx. In that case, y = Il'^{x). 

When K, = Kmax, and C : IR.,nax ^ ^max is 9- kernel operator, 11*-^ = C'* o C has 
an interpretation similar to 



where — imC* is thought of as a IR,ui,i-subscmimodule of M,uin, so that, 
P_iniC*(v) = mm{w G — imC* | w > v} . 
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where < denotes here the usual order on R . since the natural order of Mmin the 

— U — X 

reverse of the usual order. When B : M,„ax ~^ I^max is also a kernel operator, we 
have 

This factorization is instrumental in the geometrical interpretation of the finite 
element algorithm, see |AGL04I Example 10]. 

3. The max-plus finite element method 

In this section we describe the max-plus finite element method to solve Prob- 
lem Let 5* and v* be defined as in the introduction. Since 5"*+* = 5* o S** , for 
t, t' > 0, we obtain the recursive equation: 

(6) v*+^ = S^v\ t = 0,5,--- ,T-6 

with v° = (/) and S = ^ , for some positive integer N. Let W be a Rmax-semimodule 
of functions from X to R^ax such that (f) € W and for all v € W, t > 0, S*v € W. 
We suppose given a "dual" semimodule Z of "test functions" from X to Kmax- The 
max-plus scalar product is defined by {u \ v) ~ \/ xex u{x) <E) v{x), for all functions 
-u, w : X — > Mmax- We replace © by: 

(7) (z I = (z I 5^*), VzeZ , 

for t = 0, (5, . . . , T — (5, with u*, . . . , ij-^ e W. This equation can be seen as the 
analogue of a variational or weak formulation. Kolokoltsov and Maslov used this 
formulation in |KM88| to define a notion of generalized solution of Hamilton- Jacobi 
equations. We consider now a semimodule Wh C W generated by the family 
{wi}i<i<p. We call finite elements the functions Wi. We approximate -u* by e 
Wh, that is, D* ~ i;^ = 0f^;^ WiA- , where A* € Mmax- We also consider a semimodule 
Zh C Z generated by the family {zj}i<j<q. The functions zi, • • • , Zq will act as 
test functions. We replace O by 

(8) {z,\v'+') = {z,\S%), yi<j<q , 

for i = 0, (5, • • • , T - (5, with ^ cj),, ~ (j) and G , i = 0, (5, • • • , T. 

Since Equation (jS)) need not have a solution, we look for the maximal subsolution, 
i.e. the maximal solution v^^^ G Wh of 

(9a) (z, I vi+') < (z, I S'vi) yi<j<q . 

We also take for the approximate value function tjJJ at time the maximal solution 
G Wh of 

(9b) v1<v'^ . 

Let us denote by Wh the max-plus linear operator from M^^g^^ to W with matrix 
Wh = col(wi)i<i<p, and by the max-plus linear operator from W to R^i^^ whose 
transposed matrix is Zh = col(zj)i<j<q. This means that WhX = WjAj for all 

A = (Ai)i=i,...,p G K^iaxi and {^h^)j = {^j I for all w G W and j = 1, . . . ,q. 

Proposition 2 f |AGL04] ) . The maximal solution v'j^^ G Wh of ^ is given by 
'^'V'^ = Sti, where 

Si = o S' . 
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The following proposition provides a recursive equation verified by the vector of 
coordinates A* and is proved in |AGL04| . 

Proposition 3 ( |AGL04] ) . Let vj^ G Wh be the maximal solution of for t = 
0, (5, . . . , T. Then, for every t = 0^ d, . . . ,T , there exists A* € I^max such that vj^ = 
VF/iA*. Moreover, the maximal A* satisfying these conditions verifies the recursive 
equation 

(10) A* = iZlW,,)\{ZlS'W,A'-') , 
with the initial condition = Wh\(t>- 

For I < i < p and 1 < j < q, we define: 

(11) i^hh = {zj I 

(12) {Bh)r. = {z, I S'w,) 

Ah and Bh are respectively the matrices of the max-plus linear operators Z'^Wh 
and ZlS^Wh. Equation ^ may be written explicitly, for 1 < i < p, as 

AH nhn (-(A),,+ max ((i?„),fe + A*-^)) . 

This recursion may be interpreted as the dynamic programming equation of a de- 
terministic zero-sum two players game, with finite action and state spaces. 
The ideal max-plus finite clement method can be summarized as follows: 

(1) Choose 5 ^ ^ and the finite elements (u'i)i<,;<p and {zj)i<j<q, 

(2) Compute the matrix Ah by Hll|l and the matrix Bh by (|12|l . 

(3) Compute A° = Wh\(t) and = Wh\^ ■ 

(4) For t = 5,25,...,T, compute A* = A,,\(B,,A*-*) and = WhX^ ■ 

Then, approximates the value function at time t, w*. 

Fleming and McEneaney proposed a max-plus based method |FM00j , which also 
uses the linear formulation JHJ. They approximated the evolution semigroup by 
a max-plus linear semigroup acting on a finitely generated semimodulc of functions. 
A comparison of this method with the ideal max-plus finite clement method appears 
in ]ki\U)A\ . 

4. Small time approximation of the Lax-Oleinik semigroup 

Computing Ah from pi(l is an optimization problem, whose objective function 
is concave for natural choices of finite elements and test functions (see Section 
below). Evaluating every scalar product (z | S^w) leads to a new optimal control 
problem, which is simpler to approximate than Problem (Q, because the horizon 
is small, and the functions z and w have a regularizing effect. In |AGL04 | . we 
proposed to use the following approximation of S^w provided by the Hamilton- 
Jacobi equation (|2a|l : 

(13) S^wix) ~ w(x) + 5H(x, -^), for all x £ X. 

ox 

In this paper, we use the approximation of S^w by the function [S'''i(;]~ such that, 
for all x e X 

(14) [S"'u;]~(a;) = snp {siix,u) + w{x + 5 f{x,u))]. 
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Let [S^WhY denotes the max-plus linear operator from R^J^ax to W with matrix 
[S^WhY = co\{[S^WiY )i<i<p- The above approximation of S^w yields an approxi- 
mation of the matrix Bf^ by the matrix := Zj^[S^WiiY i whose entries are given, 
for 1 < i < p and 1 < j < <?, by: 

(15) {Bh)ji^ sup {zj{x) +Wi{x + 5f{x,u)) + 5l{x,u)}. 

xGX,u£U 

The following proposition shows that under assumptions on £, /, Zj and Wi, com- 
puting the approximation p5(l is a concave maximization problem. In this case, 
one can compute the entries of the matrix B^ using standard convex optimization 
algorithms. 

Proposition 4. Let X be a convex set of M." and let U be a convex set o/M'". 
Assume that z : M" R U {-oo} is concave, I G C^iX x ;7,R), w e C2(R",R) 
and f is affine. Let F = ^ , G ~ ^ and let \\ ■ \\ denotes the Euclidean norm of 
operators. Assume that there exist a,/3,C > such that —CLn < < —(iln, 
< —aim, \\^)^\\ ^ C and \ \ q^q^ \ \ < C, where /„ is the n x n identity matrix. 
Then there exists a constant Sq = So{a, P,C,\\F\\,\\G\\) .such that, for all S < 6o, 
the function {x, u) > z{x) + w(^x + Sf{x, u)) 4- Sl{x, u) is concave. 

Proof. Since z is concave, it suffices to prove that the function {x,u) t-^ b{x,u) = 
w(^x + 5f{x, u)) -I- 5i{x, u) is concave. Since w, i, f are and X and U are convex 
sets, we must show that h*V'^b{x, u)h < for all h e R"+™, x G X, u e U, where 

W^b{x,u) denotes the Hessian of b at (x,u). Set x' = x + Sf{x,u) and h 



where hi e R" and /is G M". Then 



h, 



h*W%{x,u)h = h;[6^{x) + 6^G*^{x')g') h2+ 



h 

< 



((5C(1 + 211^^11) - /?) \\hi\\^ - Sa\\h2f + 



25C{l + \\G\\+S\\F\\\\G\\)\\hMh: 



2 



We recognize a quadratic form in the variables and ||/i2||- This quadratic 

form will keep a constant sign (which is negative) for non zero values of h if its 
discriminant is negative, i.e. if 

S^C^\\Ff\\Gf + 26^C^\\F\\\\G\\{1 + \\G\\) + 

6C(^C{1 + \\G\\f + a{l + 2||F||)) < a/3. 

This is the case in particular if 

3maxf<53c2||F|n|Gf,252^2|j^|j|jg||(^^||g||)^ 



SG 



(c(l + ||G||)' + a(l + 2j|i^||))] <aP, 
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and so, we can take 

af3 / a/J 



So = min ( i 



3C2||F|12||G|i2' Y 6C2||F||||G||(l + ||G| 
af3 



3c[c{l + \\G\\j' + a{l + 2\\F\\)) 

□ 

5. Error analysis 

We first recall a general lemma showing that the error of the finite element 
method is controlled by the projection errors, ||n^^^z;* — v*\\oo and ||n'^'>i;* — w*|ioo, 
and by the approximation error, — S'^Wiljoo- 

Lemma 5 f )AGL04j ). Fort = 0, (5, • • • ,T, let be the value function at time t, and 
t;*j be its approximation given by the max-plus finite element method, implemented 
with the approximation B^^ of B^, given by I^15\) . We have 

IK',T-^'^lloo < (1 + T)f max \\[S^w,T - S^w,\\^ 

\ l<i<p 

+ sup (||n^'*«*-i>loc + ||nv^^v*-i;*|ioo) 

t=0,A\...,T 

The proof of this lemma uses the fact that projectors over max-plus semimodules 
are non-expansive in the sup-norm. 

To state an error estimate, we fix a norm || • || on K." and wc make the following 
assumptions: 

{HI) The semigroup preserves the set of i-semiconvex functions, for some c > 0. 
{H2) f : X X U — > R" is bounded and Lipschitz continuous with respect to x: 
3i/ > 0, Mf > such that 

\\f{x,u)- f{y,u)\\<Lf\\x-y\\ yx,yex,yueu 

\\f{x,u)\\ < Mf, yxeX,ueU. 

{H3) £ : X X U — > M is bounded and Lipschitz continuous with respect to x: 
3Li > 0, Ml > such that 

\e{x,u)-e{y,u)\ <Le\\x-y\\ \fx,yeX,ue U, 

\e{x,u)\<Mi, yxeX,ueU. 

(HA) (f> : X ^M. is bounded and Lipschitz continuous. 
Recall that a function / is ^-semiconvex if ^HxH^, where || • II2 is the standard 

Euclidean norm of R", is convex. Spaces of semiconvex functions were already used 
by Fleming and McEneaney jFMOOj . 

Wc shall use the following finite elements. 

Definition 6 (Lipschitz finite elements). We call Lipschitz finite element centered 
at point X £ X, with constant A > Q, the function w{x) = —A\\x — x\\i where 
\\x\\i ~ X]r=i l^-norm o/M". 

The family of Lipschitz finite elements of constant A generates, in the max-plus 
sense, the semimodule of Lipschitz continuous functions of Lipschitz constant A 
with respect to |j • ||i. 
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Definition 7 (Quadratic finite elements). We call quadratic finite element centered 
at point X d X , with Hessian > 0, the function w{x) = —-^Wx — x\\2. 

When X ~ M", the family of quadratic finite elements with Hessian i gener- 
ates, in the max-plus sense, the scmimodule of lower-semicontinuous i-scmiconvex 
functions. 

Lemma 8. Let X be a compact o/R". We make assumptions (H2) and (H3). 
Assume that w and its derivative are both Lipschitz continuous. Then there exists 
Ki>0 such that WiS^wf - S^w\\oo < KiS^- 

Proof. Denote by Mdw and Md2^ the Lipschitz constants with respect to norm 
II • II of ui and its derivative respectively. Wc first show that there exists K > such 
that S^w{x) - [S^wf (x) > -KS'^. For all x € X we have 

iS^w){x) > sup I /p £ix{s), u)ds + w{x{S)) \u e U, x{0) = x, 

x{s) = f{x{s),u),se [0,5]}. 

In other words, we bounded S*w from below by considering only constant controls. 
Since 

[^(a;(s), m) — ^(x, u)](is < L^ \\x{s) — x\\ds 

Jo 

.5 



< Li [ sMfds 
Jo 



1 
2 

we obtain 



< ;:Lf MfS^ 



X2 , 

{S^w){x) > -LiMf— + sup |5€(.T, u) + w{x{6))\u € U, x{Q) = x, 

x{s) = f{x{s),u),se [0,(5]}. 

Moreover 

\w{x{S)) - w{x + 6f{x,u))\ < MdwW / f{x{s),u)ds-df{x,u)\\ 



< MDn,LfMf- . 

We deduce that 

^2 

{S^w){x)-[S^wr{x) > -{Li + MD^Lf)Mf— . 
We now prove an opposite inequality. For x Cz X we have 

£{x{s),u{s))ds + w{^x + / f{x{s),u{s))ds^ \ 

u{s) e U,x{s) = f{x{s),u{s)),x{0) ^ x,se [0,(5]}. 
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By the same arguments as before, we show that 



^2 

{S^w){x) < {Le + MDwLf)Mf—+ 



sup 

M(s)ec/,Vs Jo 



£{x,u{s))ds + w(^x + / /(x, 'ii(s))rfs) I 



Using the fact that 

w{x + J f{x, u{s))ds) < w{x) + Vw{x) ■ J f{x, u{s))ds^ + Af£,2„M^ 

w{x + Sf{x, U)) > W{x) + Vwix) ■ (Sf{x, m) j - Mo^njMf 



5_2 
2 

and 

[e{x, u{s)) + ^w{x) ■ f{x, u{s))] ds < sup Si{x, u) + Vw(a;) • {Sf{x, u)) , 

u£U 

we deduce that 

S^w{x) - [S^wf {x) < {Le + MDu,Lf + 2MD--n,Mf)Mf — . 

□ 

Using Lemma |5| Lemma [S] and exphcit estimates of the projection errors ap- 
pearing in Lemma|5l along the hnes of |Lak03j . we derive the fohowing convergence 
result. 

Theorem 9. Let Q, he an open convex set ofW^, X = VL and X = X ^ B(0, ci), 
where L > 0, c > 0. Suppose that X and X have regular grids T/^x o,nd Taz: 
respectively of size Ax. We make assumptions (H1)-(H4), and assume that the 
value function at time t, w*, is L-Lipschitz continuous with respect to \\ ■ ||i and 
^-semiconvex for all t > 0, with the same constant c as in (HI). Let us choose 
quadratic finite elements Wi of Hessian ^ , centered at the points of Tak . Let us 
choose, as test functions Zj, the Lipschitz finite elements with constant A > L, 
centered at the points of T/^^. For t = 0, 5, . . . ,T , let w*^ be the approximation of t;* 
given by the max-plus finite element method implemented with the approximation 
B'^ of Bh given by H15|) . Then, there exists a constant K > such that 

Remark 10. A different approximation of Bh relying on was used in |AGL04| 
ILak03 | . It is easier to implement. In particular, it avoids the numerical solution of 
optimization problems, as soon as explicit formulcE are available for the Hamiltonian 
H and the point of maximum of {x) + Wi{x), which is frequently the case. 

However, it only leads to an error of order \/S + ^ . 
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6. Numerical results 

Wc now present some results obtained by the method discussed above. 

Example 11 (Linear Quadratic Problem). We consider the case where [/ = R^, 
X = R'^, (j>=0, 

£{x, u) ^ '^^'^ _ Uj^+U2 f{x,u)^u . 

We choose quadratic finite elements Wi and Zj of Hessian -. We represent in 
Figure^ the solution given by our algorithm in the case where T — 5, S — 0.5, 
Ax = 0.05, c = 0.1. 



Figure 1. Max-plus approximation of a linear quadratic control 
problem fExample lll() 

We observe a truncation effect on the boundaries of the set X . If we restrict X 
to the set [—0.5,0.5]^, we obtain a Loo-error of order 0.07. 

Example 12 (Distance problem). We consider the case where T = 1, = 0, 

X=[-l,l]2, C/=[-l,l]2, 

£{x, u) = 




f{x,u) = 

We choose quadratic finite elements Wi of Hessian ^ and Lipschitz finite elements 
Zj with constant A. We represent in Figure\^the solution given by our algorithm 
in the case where T — 1, S — 0.1, Aa; = 0.1, A = 3 and c = 1. The Loo-error is of 
order 0.15. 
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